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A First Look at the Auriga- California Giant Molecular Cloud 
With Hers cheF and the CSO: Census of the Young Stellar Objects 

and the Dense Gas 

Paul M. Harvey 1 , Cassandra Fallscheer 2 , Adam Ginsburg 3 , Susan Terebey 4 , Philippe 
Andre 5 , Tyler L. Bourke 6 , James Di Francesco 7 , Vera Konyves 5 ' 8 , Brenda C. Matthews 7 , & 

Dawn E. Peterson 9 

ABSTRACT 

We have mapped the Auriga/California molecular cloud with the Herschel 
PACS and SPIRE cameras and the Bolocam 1.1 mm camera on the Caltech 
Submillimeter Observatory (CSO) with the eventual goal of quantifying the star 
formation and cloud structure in this Giant Molecular Cloud (GMC) that is 
comparable in size and mass to the Orion GMC, but which appears to be forming 
far fewer stars. We have tabulated 60 compact 70/160 /im sources that are 
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likely pre-main-sequence objects and correlated those with Spitzer and WISE 
mid-IR sources. At 1.1 mm we find 18 cold, compact sources and discuss their 
properties. The most important result from this part of our study is that we find 
a modest number of additional compact young objects beyond those identified at 
shorter wavelengths with Spitzer. We also describe the dust column density and 
temperature structure derived from our photometric maps. The column density 
peaks at a few x 10 22 cm -2 (iV# 2 ) an d is distributed in a clear filamentary 
structure along which nearly all the pre-main-sequence objects are found. We 
compare the YSO surface density to the gas column density and find a strong 
non-linear correlation between them. The dust temperature in the densest parts 
of the filaments drops to ~ 10K from values ~ 14-15K in the low density parts of 
the cloud. We also derive the cumulative mass fraction and probability density 
function of material in the cloud which we compare with similar data on other 
star-forming clouds. 

Subject headings: infrared: ISM — stars: formation — ISM: clouds — ISM: 
structure — ISM: individual objects: Auriga- California GMC 



Introduction 



The Auriga-California molecular cloud (AMC) is a large region of relatively modest 
star formation that is part of the Gould Belt. We have adopted the name 'Auriga- California 



Molecular Cloud" since the region is listed as 'Auriga" in the Spitzer Space Telescope (Werner 



et al. 2004) Legacy Survey by L. Allen, while it has been called the "California Molecular 



Cloud" by Lada, Lombardi & Alves (2009) based on its proximity to the "California Nebula". 



The Spitzer observations of this region are described by H. Broekhoven-Fiene et al. (2013, in 
preparation) as part of the large scale Spitzer "From Cores to Planet-Forming Disks" (c2d) 
and "Gould Belt" programs that were aimed at cataloguing the star formation in the solar 
neighborhood. A similar large-scale mapping program with the Herschel Space Observatory 



(Pilbratt et al. 2010), the "Herschel Gould Belt Survey" (KPGTl_pandre_01) (Andre et al. 



2010), has been observing most of the same star-forming regions, but the AMC was not 



included in the original target list for that program. 



* Herschel is an ESA space observatory with science instruments provided by European-led Principal 
Investigator consortia and with important participation from NASA. 
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The AMC provides an important counterpoint to other star-forming regions in the Gould 



Belt, particularly the well-known Orion Molecular Cloud (OMC). As described first by Lada, 



Lombardi & Alves (2009), the AMC is at a likely distance of 450 pc (though Wolk et al. 



(2010) quote a slightly larger distance of 510 pc). This distance is quite comparable to that 



of the OMC, and the mass of the AMC estimated by Lada, Lombardi & Alves (2009) is also 
quite similar, ~ 10 5 M . The most massive star that is forming in the AMC, however, is 
probably the Herbig emission-line star LkHalOl, likely an early B star embedded in a cluster 



of lower mass young stars (Andrews & Wolk 2008 Herbig et al. 2004). This situation is in 



stark contrast to the substantial number of OB stars found in several tight groupings in the 
OMC ( Blaauw|1964 ). Lada, Lombardi & Alves (2009) investigated the distribution of optical 
extinction in the AMC and used those results together with 12 CO maps from Dame et al. 



(2001) to conclude that one possibly significant difference between the AMC and OMC is 



the much smaller total area exhibiting high optical extinction in the AMC, roughly a factor 
of 6 smaller area above Ak = 1 mag. 

Herschel observations have demonstrated probably the best combination of sensitivity 
and angular resolution to a range of dust column densities in star-forming regions, as well 
as excellent sensitivity to the presence of star formation from the very earliest stages to the 
so-called Class II objects with modest circumstellar disks. We therefore have undertaken 
a Herschel imaging survey of a ~ 15 deg 2 area of the AMC to document the full range in 
evolutionary status of the star formation in this cloud as well as the distribution and column 
density of dust as a proxy for the total mass density. We have supplemented the Herschel 
observations with a 1.1 mm Bolocam map from the Caltech Sub millimeter Observatory 
(CSO) to identify the extremes in cold, dense material. We describe the observations and 
data reduction in the following section. Then in §|3] we discuss our extraction of the compact 
source component in the 70/160 fim Herschel data as well as in the 1.1 mm maps and 
compare our fluxes with those from other measurements. In ^4]we describe several interesting 
individual objects. In §[5] we discuss the dust column density and temperature maps derived 
from our Herschel PACS/SPIRE images and the relationship between this dust emission and 
previous observations of dust absorption and gas emission. We also derive a quantitative 
correlation between the gas density and YSO surface density. Finally in ^6] we begin a 
discussion of the differences between star formation in the AMC versus that in the OMC, a 
subject which we will investigate more fully in future studies. 
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Observations and Data Reduction 



2.1. Herschel Observations 



Our Herschel program, the "Auriga- California Molecular Cloud" (OTl_pharve01_3), was 
designed to use the same observing modes as comparable parts of the large-scale Gould Belt 



program by Andre et al. (2010). For both programs the "Parallel Mode" of PACS/SPIRE 



(Griffin et al. 2010 ) was used to cover the largest possible size region in a reasonable observing 



time, and a much smaller region was covered with PACS (Poglitsch et al. 2010) alone to 



provide additional sensitivity and wavelength coverage. The Parallel Mode observations 
were done with PACS at 70 jum and 160 /im, and the SPIRE observations naturally included 
the three SPIRE photometric bands, 250 /im, 350 /im, and 500 /im, that are observed 
simultaneously. With the PACS-only observations, as for the larger scale Herschel Gould 
Belt program, we used PACS at 100/160 /im with slow scan speed (20"/s) which essentially 
preserves the full diffraction-limited resolution of Herschel. These latter observations were 



centered on the well-known LkHalOl cluster (Andrews & Wolk 2008) which includes a 



significant fraction of all the obvious star formation in this cloud. We do not discuss these 
PACS-only observations further in this paper, but they will be used in a subsequent study to 
help address source confusion in the dense central cluser. The total area covered in Parallel 
Mode is 18.5 deg 2 with 14.5 deg 2 covered with overlapping perpendicular scans for good drift 
cancellation. Figure [T] shows the area covered in Parallel Mode overlaid on the extinction 



map of a much larger portion of this area discussed by Dobashi et al. (2005). Our covered 



area was chosen to include essentially all of the high-extinction parts of the cloud with the 
exception of L1441 which is beyond the right (low Galactic longitude) end of our maps. The 
PACS-only observations covered 1.4 deg 2 . The details of the observations and ObsIDs are 
listed in Table [1} The observed Parallel-Mode area was divided into three separate pieces 
for efficiency in AOR design and observatory scheduling. The area covered includes nearly 
all of that observed by the Spitzer Gould Belt study of the AMC with the exception of a 
small separate portion northwest of the end of our maps. 

The initial data reduction process is essentially the same as that used for several other 



star-forming clouds from the Gould Belt Survey, e.g. (Sadavoy et al. 2012 Peretto et al. 



2012). The first step consists of reducing the Herschel data to level 1 products using the 



Herschel Interactive Processing Environment (HIPE) version 8.1.0 (Ott 2010). Maps of the 



three sub- regions listed in Table [T] were obtained using Scanamorphos version 16 (Roussel 



2012) using the two perpendicular scanmaps to remove correlated noise such as low frequency 



drifts. The pixel scales for these maps were 3.2", 5", 6", 10", and 14" respectively at 70 /un, 
160 /zm, 250 /im, 350 /im, and 500 /im. These individual maps are shown in Figures |2]-[6 
(electronic edition only). We then used two different source extractor routines. The first 
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Fig. 1. — Extinction image of Dobashi et al. (2005) with the outline of our covered area 



shown to illustrate that we have observed most of the high extinction parts of this cloud. 
The maximum extinction is = 4.5 mag and the minimum is essentially zero. The image 
is oriented in Galactic coordinates and covers an area of 12° (/) x 7.3° (b). 
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was the getsources package (version 1.120526) QMen'shchikov et al.|2012| ) that was developed 
to search for sources over a range of spatial scales and extracts sources simultaneously over 
multiple bands that have substantial differences in angular resolution. The second source 
extractor was the c2dphot package developed as part of the Spitzer Legacy c2d program 
( |Harvey et aL]|2006[ |Evans et al. | |2007[ ) which was designed to work with point-like and small 
extended sources up to roughly twice the beam size and was based on the earlier DOPHOT 



package (Schechter, Mateo, & Saha 1993). In this paper, we make use mostly of the results 



from the c2dphot processing (shown in Table [2]) since we are primarily addressing point-like 
and very compact sources (§[3j in addition to the very large scale structure (§{5]). Future 
publications will use the results of the getsources processing to investigate the medium-scale 
emission. 

Figure [7] shows a 3-color composite (70 /im, 160 /im, and 250 jttm) of the entire region 
mapped at 70 /im, 160 /im, 250 /im, 350 /im, and 500 /im. The two most obvious features 
of this map are: 1) the bright collection of sources and nebulosity at the left end of the map 
(Southeast) where the LkHalOl cluster is located, and 2) the long network of filamentary 
structure that pervades much of the mapped area. Such filamentary structure is now known 
to be typical in Galactic star-forming regions from the work of the Herschel Gould Belt 



survey (Andre et al. 2010) as well as the Herschel Galactic Plane Survey, HIGAL (Molinari 



et al. 2010) and has also been discussed earlier by Myers (2009). Subsections of some of the 



mapped areas are discussed in more detail later in £|4| 

In addition to the basic map-making and source extraction, we also present in Figure 
[8] results on dust temperature and optical depth over the entire mapped area. We used a 



method similar to that described by Konyves et al. (2010); we first determined zero-point 



offsets following the procedure described by Bernard et al. (2010) and then convolved the 



shorter wavelength Herschel images to the resolution of the 500 /im data. We derived SED 
fits to the 160 /im, 250 /im, 350 /im, and 500 /im data for each pixel position in the maps 
using a simplified model of dust emission, F u — k v x B(u,T)x column density. We assumed 
a dust opacity law of k v = 0.1(b / /1000GHz) l3 cm 2 / g) and fixed the dust emissivity index 
to (5=2 with the standard mean molecular weight, /i = 2.33. Because of the high S/N 
over most of the area of the flux maps that were used to derive these column-density and 
temperature maps, the major uncertainty in the absolute values of T and r are those due 
to the inherent assumptions in using the equations above to represent the dust emission. It 
is likely, though, that the maps provide an excellent representation of relative temperatures 
and column densities with absolute uncertainties of order ± 15-20% in temperature and ± 
a factor of 2 or more in column density. We discuss these maps more fully in ^5] where we 
compare the column densities to those derived from extinction measurements and analyze 
the distribution of star formation relative to the inferred gas densities (note, all column 




Fig. 3.— False color image of the 160 /im map derived with Scanamorphos from our 
Parallel-Mode observations. 



-8- 




Fig. 5.— False color image of the 350 /im map derived with Scanamorphos from our 
Parallel-Mode observations. 
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Fig. 6.— False color image of the 500 fim map derived with Scanamorphos from our 
Parallel-Mode observations. 
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Fig. 7. — False color image with 70 /iiri (blue), 160 /zm (green), and 250 jum (red) of the 
mapped area. The locations of the sources in Tableware marked with a red 'X', and those 
that are not in the Spitzer YSO list of H. Broekhoven-Fiene et al. (2013, in preparation) are 
also surrounded with a black square. 
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densities discussed in this paper are measured as N H2 ). 



The digital versions of all these maps will be available soon after publication of this 
paper on the Herschel Science Center's web site for user-provided data, 



http:/ /herschel. esac.esa.int/UserProvidedDataProducts.shtml 



2.2. CSO Observations 

We used the Bolocam imaged at a wavelength of 1.1 mm to map much of the area 
covered in our Herschel observations during the nights of 14-16 November 2011. We utilized 
observing techniques similar to those used for the Bolocam Galactic Plane Survey (BGPS) 



as described by Aguirre et al. (2010) and Ginsburg et al. (in prep.). Alternating maps were 
made scanning roughly parallel and perpendicular to the Galactic plane at a scan speed 
of 120"s _1 . Multiple overlapping maps were obtained over roughly the eastern 2/3 of the 
Herschel mapped area; the total area observed was 6 deg 2 ; the area covered is indicated 
in Figure [8j Due to non-uniform coverage and varying weather conditions the noise in the 
Bolocam maps is not constant, but is typically ~ 0.07 Jy/beam. This is substantially higher 



than the noise in maps of several other Gould Belt clouds presented by Enoch et al. (2007) 



0.01-0.03 Jy/beam, due to our significantly smaller observing time per pixel. The primary 
flux calibrator was Uranus. The map data were reduced using the software described by 



Aguirre et al. (2010) for the Bolocam Galactic Plane Survey (BGPS), utilizing correlated 
sky-noise reduction with 3 PCA (Principal Component Analysis) components. Following 



that, sources were extracted as described by Rosolowsky et al. (2010) for the BGPS. In 



addition to this large-scale mapping, we also observed a small area centered on one of the 
strong Spitzer sources to the northwest of the scanned region, SSTGB04012455+4101490, 
for which we have no corresponding Herschel data. 



Aguirre et al. (2010) have carefully investigated the inherent spatial filtering that occurs 



in removing correlated sky-noise in ground-based observations at this wavelength. For the 
case of subtraction of 3 PCA components roughly half the flux is lost for structure larger 
than 300". Indeed, the largest coherent area of 1.1 mm emission in our map is a 4' wide 
area centered on LkHalOl. Therefore, we present the results from this part of the study as 
positions and flux densities for the compact emission regions detected. Table [3] lists these 
positions and the fluxes within several different apertures for the compact sources detected 
at 1.1 mm with peak S/N > 2. Note, Source 2 is the bright galaxy 3C111 which was also 



^http://www.cso.caltech.cdu/bolocam 
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Fig. 8. — Upper image: column density image with positions of Bolocam 1.1 mm sources 
marked with squares (3C111 is indicated with the diamond). The area covered by our 1.1 mm 
Bolocam mapping is outlined with the dash-dot line. The highest column areas (white) have 
N_ff2 ~ 5 x 10 22 cm -2 while the lowest column areas represent values ~ 1 x 10 21 cm' 2 . Lower 
image: dust temperature image with maximum T^ ~ 28K at LkHalOl near the left end of 
the image and minimum temperatures of order 10K in the darkest parts of the filaments. 
The median derived dust temperature over most of the area is ~ 14.5 K. 
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our primary pointing calibrator and secondary flux calibrator. Table [3] also lists the Herschel 
sources from Table [2] that are located within 45" of each 1.1 mm source position and likely 



associated with it. We discuss these 1.1 mm sources below in £3.2 



3. Compact Sources 

3.1. The 70 /im Objects 

The goal of our investigation of the compact sources in the AMC is to complete the search 
for pre-main-sequence and proto-stars that began with the Spitzer Gould Belt program (H. 
Broekhoven-Fiene et al. 2013, in preparation) and, in particular, to search for the most dust- 
enshrouded objects that might have been missed by that program because they emit most of 
their luminosity in the far-IR. Herschel at 70 /im provides the highest resolution imaging in 
the far-IR of any current or planned facility, and conveniently the 70 jum resolution (X/D ~ 
4") is also nearly identical to that of Spitzer at 24 /im. Although the resolution of Parallel- 
Mode observations is not quite as high as Herschel's diffraction limit because of image blur 
from the fast scan speed in Parallel-Mode, the resolution achieved is not much below that 
limit. Therefore, an additional goal of this investigation is to use this resolution to measure 
fluxes in the far-infrared more reliably than Spitzer in confused regions. With a complete 
and reliable census of all the stages of star-formation in the AMC, we will be able to make 
the most informative comparison of it with the OMC 

The source extractor c2dphot operates in two modes. In the first mode, it searches 
through the image at sequentially lower flux levels for local maxima, characterizes them 
as point-like or extended (ellipsoidal), and subtracts them from the image. In the second 
mode the code is given a list of fixed positions at which it fits the point-source-function 
(PSF) to whatever flux above the background exists at that position. This mode is useful 
for determining upper limits and for testing for faint objects in the wings of bright ones. In 
both modes an aperture flux is calculated as well as the PSF- or ellipsoidal shape-derived 
flux. To find the most complete set of possible objects to correlate with objects at other 
wavelengths, we first processed both the 70 /im and 160 /im images in the most general 
c2dphot mode, allowing the code to fit flux, position, and shape down to the lowest flux 
levels present in the image, i.e., essentially the noise level. This process produced a list of 
~6500 sources at 70 /im and ~500 sources at 160 /im, of which probably over half are noise at 
both wavelengths. After comparing a number of individual cases while trying to correlate the 
160 /im objects with those found at 70 /im, we identified two complicating issues. First, the 
obvious issue of the larger PSF at the longer wavelength meant that sometimes more than one 
70 /xm source would be within the 160 /im PSF. Second, because cooler, more extended dust 
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is naturally detected at the longer wavelength, in some cases the 160 /im source equivalent to 
a nearly point-like 70 /im source would be extended and have an asymmetric shape, making 
an automated detection and association difficult. For these reasons we decided to determine 
the 160 /im fluxes (or limits in most cases) for the 70 /im detections by running c2dphot in its 
second, fixed-position mode at 160 /im, using the 70 /im detection list for the input positions. 
In this case it can also be useful to compare the PSF-fit fluxes with those determined from 
aperture photometry as a secondary indication of larger source extent or confusion. For the 
six coldest sources discussed later in Figure [14] with < 40K we have also extracted flux 



densities at 250 /im in the same way as the 160 /im fluxes and show the PSF-fit fluxes. Since 
the 250 /im PSF has a full-width-half-maximum of roughly 18", we have not measured the 
fluxes of the bulk of our objects beyond 160 /im; the issues with assigning fluxes to individual 
sources at 160 /im are incrementally more problematic at the SPIRE wavelengths. A future 
study making use of the getsources processing is likely to produce the most reliable long 
wavelength SEDs for most of the sources. 

The nominal absolute calibration uncertainty in PACS photometry is now believed to be 
±3 % at 70/im and ± 5% at 160/iirj^J These values only apply to well-sampled, color-corrected 
point sources. The SPIRE photometry is believed to be calibrated to ± 5% under equally 
ideal circumstances^] For the purposes of this study, we assume the more conservative value 



of ± 15% as used, for example, by Konyves et al. (2010) from the original instrument papers 



by Poglitsch et al. (2010) and Griffin et al. (2010). Given that we have flux measurements 



for a number of sources at 70 /im also from the Spitzer Gould Belt program as well as 
flux determinations from our Herschel data set using the completely different getsources 
algorithm, we have an opportunity to check our flux measurements for systematic effects 
and problems like non-linearity. Figure [9] shows plots of these two comparisons. For the 
Spitzer comparison we used the list of reliable YSOs described by H. Broekhoven-Fiene 
et al. (2013, in preparation) and excluded several objects in very confused regions. The 



Spitzer photometry was produced by the version of c2dphot described in detail by Evans 



et al. (2007) for the final delivery of c2d data to the Spitzer Science Center. For the best 



getsources comparison (red diamonds in Figure [9]) we used only sources found to be isolated, 
well-fit by a point source at 70 /im, with a product of major axis and minor axis less than 
150" 2 , a total flux less than 1.5 x the point-source flux, and a S/N > 10. The mean ratio 
of Herschel to Spitzer fluxes is 1.03±0.3, and the mean ratio of c2dphot fluxes to getsources 
fluxes is 0.96±0.13 for the sources marked with red diamonds. The scatter between the 
Herschel and Spitzer results appears to be independent of flux level, while that between 



^http://herschel.esac.esa.int/twiki/bin/view/Public/PacsCalibrationWeb?template=viewprint 
§ http://herschel.esac. esa.int/hcss-doc-9.0/ 
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the two methods used on the Herschel data is consistent with what one would expect as 
a function of the signal-to-noise level. This excellent agreement, particularly for the two 
methods used on the Herschel data, between flux determinations over a very wide range of 
brightness suggests that both c2dphot and getsources provide highly reliable extractions and 
flux determinations, certainly for compact objects. 

To identify the objects in the AM C that are most likely to represent young sources with a 
stellar or pre-main-sequence core, we culled our starting list of ~ 6500 sources in several ways 
using shorter wavelength data. Much of our observed area has been covered at 24 /xm in the 
Spitzer Gould Belt survey of H. Broekhoven-Fiene et al. (2013, in preparation); in areas not 
observed with Spitzer the recently released WISE all sky survey ( Wright et al.||2010 ) provides 



relatively deep mid-infrared photometry over a similar range of wavelengths. To start our 
search for reliable young objects we included only those sources with S/N > 7.5 and F v at 70 
/im> 85 mJy with at least a 2.5<r detection at one of the closest neighboring wavelengths, i.e. 
22/24 /iin Spitzer-MIPS/WISE, or our 160 fim Herschel photometry. These criteria reduced 



the list of possible young objects to 513 sources. Koenig et al. (2012) have estimated a 
contamination rate for extragalactic sources in star-forming regions at the sensitivity and 
wavelength of the WISE survey of ~ 10 objects per square degree. This high background 
level suggests that a significant fraction of our 513 candidate sources is extragalactic. 

To eliminate as many extragalactic sources from the candidate list as possible, we used 
a multi-pronged approach that relied on: (1) examination of individual images for likely 
galaxies at 70 jttm as well as 24 /zm (Spitzer), 2MASS, and the red (DSS) Digitized Sky 



Survey images, (2) Spitzer /WISE color-color diagrams using the criteria developed by Koenig 



et al. (2012), (3) the 2MASS "gaLcontam" flag which signifies a likely extended extragalactic 
object (gaLcontam = 1), and (4) objects of low S/N in the images or those that had no clear 
point-like core at 70 /xm. Because of the wide range in colors, brightnesses, and angular sizes 
of the extragalactic objects, all of these criteria contributed substantially to the elimination 
process. After this triage we were left with 209 possible young candidates. This sample was 
clearly still "polluted" with extragalactic sources and evolved stars as we found by searching 
the SIMBAD data base and examining several dozen sources in DSS images. The most 
likely YSO candidates found by Spitzer (H. Broekhoven-Fiene et al. 2013, in preparation), 
however, are located exclusively in the areas of high column density illustrated in Figure [HJ 
as are the 1.1 mm sources found in our Bolocam survey. Most of our 209 70 fim candidate 
young objects, though, are distributed quite uniformly over the area, as were the previously 
culled extragalactic objects and candidates. Therefore, at the risk of missing a likely very 
small number of young objects outside the high-column-density areas, we decided to apply 
one more criterion to our search list, i.e., to require the column density, as measured by 
Herschel at the 35" resolution of the column density map, to be above 5 xlO 21 cm -2 (Nj^) 
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Fig. 9. — Comparison of 70 /iiri fluxes (mJy) derived in this study. Left panel: ratio of the 
fluxes from the Spitzer Gould Belt survey to the Herschel c2dphot fluxes from this study for 
the YSOs in common. Right panel: ratio of the fluxes using the getsources algorithm with 
our Herschel maps to those from our c2dphot processing for the high S/N objects found not 
to be extended on the basis of the quality of fit or derived source size. Those marked with 
red diamonds also included a criterion on ratio of total flux to PSF-fit flux being less than 
1.5. 
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at the source position. This threshold reduced the young candidate list to 60 objects. To be 
confident that this criterion did not eliminate any cold, dense young sources, we examined 
the DSS and 2MASS images of the few objects in lower column density regions that had 
Fwo/F m > 1.5 and F 70 > 150 mJy, and all appeared extragalactic, i.e. not point sources. 
Therefore, the only young objects that we may have missed would be relatively faint and blue. 
Indeed, within the sample of Spitzer YSOs (H. Broekhoven-Fiene et al. 2013, in preparation) 
only one out of the ~60 objects in the Class I — II range lies outside the N# 2 > 5 x 10 21 cm -2 
area of our column density map. The combination of all of the above criteria for eliminating 
extragalactic objects was important in reaching our final sample; for example, if we had 
simply applied the column density criterion alone, we would have extracted a sample of 
~ 120 objects, half of which obviously would have been background extragalactic sources 
behind higher column density local material. As an example of the contamination issue from 



extragalactic objects, Figure 10 shows a small field in the filament north of LkHalOl with 
6 YSOs identified by H. Broekhoven-Fiene et al. (2013, in preparation) (squares) and 4 
objects identified as extragalactic from the 2MASS "gaLcontam" flag (diamonds) that are 
also extended when examined carefully at 70 /im. 



Table [2] lists the positions and 70 /xm/160 /mi c2dphot flux determinations of the final 
list of 60 objects found at 70 /mi that appear to be reliable young members of the AMC. The 
uncertainties listed are the statistical uncertainties of the measurements only; the absolute 
calibration uncertainties of ±15% have already been mentioned. For those sources that 
are identified also in the Spitzer Gould Belt data set, the 24 /im Spitzer flux is given, and 



if that is not available, the 22 /im WISE (Wright et al. 2010) flux is listed if available. 



The objects that are identified as YSO candidates by H. Broekhoven-Fiene et al. (2013, in 
preparation) are indicated, and where possible we list the object type shown by SIMBAD. 
Note, we think it unlikely that the 4 objects listed as "PN?" are truly evolved objects based 
on both their photometry and location in the AMC cloud. Table [2] also lists the spectral 



slope a (a 



dlog(\F(\)) ' 
dlog(X) 



determined from all existing members of the set of 3.6 /im, 24 



/xm, 70 /tm, and 160 /im photometry^] also shown is the corresponding YSO class based on 
the nomenclature of Lada| ( |1987[ ) as extended by Greene et al. (1994) as well as the total 



luminosity and the bolometric temperature as defined by Myers & Ladd (1993) determined 



over the same 3.6 /tm-160 /im wavelength range. We have extended the classification to 



"Class 0" to signify the most dust-enshrouded objects as suggested by Andre et al. (1993). 



Since we do not have high angular resolution photometry for most of the objects at A ~ 1 mm, 
we have used the bolometric temperature to identify these objects within the nominal Class 



^Note: the original definition of "a" used photometry only out to 24 /im, so our spectral slopes are not 
directly comparable to earlier measurements in many cases. 



Fig. 10.— Field (10' x 13') centered at J2000 R.A. 04 h 28 m 33 s , decl. +36° 25' 20" with 
6 sources marked on the 70 /im image from Table [2] with squares and their source number, 
and 4 objects with diamonds identified as extragalactic using the 2MASS "gaLcontam" flag. 
Orientation is North-up, East-left. 
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I category as defined by spectral slope. We used the criteria that objects with T bo i < 50 K 
are likely candidates for Class objects, and those with 50 K <Tt> i < 70 K we have marked 
as Class I/O, since some of them are likely to be identified as Class when the requisite 
submillimeter photometry exists and reliable submillimeter to bolometric luminosity ratios 
can be derived. The class determinations are generally similar to those found for the YSO 
candidates of H. Broekhoven-Fiene et al. (2013, in preparation), though the addition of our 
more reliable 70 /im and 160 fim photometry has made changes for a few. 

H. Broekhoven-Fiene et al. (2013, in preparation) report a total of 164 YSO candidates 
based on Spitzer and WISE photometry at A < 24/rni within the area covered by our Herschel 
survey. Clearly a significant number of these are not detected in our study. We have examined 
this list of non-detections and find that most of them are simply too faint and blue to be 
likely to be detected in the far-infrared at our sensitivity level. Several redder and brighter 
Spitzer YSO candidates lie within the confused region around LkHalOl. 

There are 11 objects in Table [2] that are not in the Spitzer YSO candidate lists of 
H. Broekhoven-Fiene et al. (2013, in preparation). Two of these (Sources 4 and 42) are 
completely undetected at 22/24 fim by WISE / Spitzer. We were able, however, to derive 
rough upper limits from the existing MIPS 24 jum images of 2 mJy for source 4 and 3 mJy 
for Source 42. Clearly both these objects exhibit relatively cold spectral energy distributions 
(SEDs). Most of the remaining 9 objects not found as YSO candidates by H. Broekhoven- 
Fiene et al. (2013, in preparation) were not selected by them simply because at least one 
IRAC or WISE band was missing from the detection list making classification as a YSO 
impossible using the Spitzer c2d/GB criteria. One of the 11 objects not identified previously 
as a YSO is Source 49 which has a very blue SED in the IRAC bands, but very strong 



far-infrared emission. It is one of only 6 objects with T bo i < 40K (see also Figure 13 and §4|. 



3.2. The Bolocam Sources 

Figure [8] shows the location of the 18 1.1 mm sources from Table [3] that are within 
our Herschel area relative to the column density derived from our Herschel mapping. With 
the exception of the galaxy 3C111 all of the 1.1 mm peaks fall on filaments of high column 
density. Excluding 3C111 and LkHalOl, 3/4 of the 1.1 mm sources have Herschel 70 fim 
objects associated with them, and these are generally from the two "earliest" YSO classes, 
and I (Note: Bolocam source 1 is associated with 70 /im source 9, which is a very bright 



emission-line star that is a likely FU Orionis object (Sandell & Aspin 1998)). There are, 
however, four 1.1 mm emission sources with no associated Herschel source from Table [2] 
Source 10 is associated with a diffuse "blob" of emission near LkHalOl in our 70 /im map, 
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but Sources 3, 4, and 5 have no clear 70 /iin counterpart. They are also three of the 1.1 
mm sources with the lowest S/N ratio, but as just noted, they are clearly located on high- 
column-density filaments. 

If the bulk of the 1.1 mm emission from the sources in Table [3] arises from very cool, 
dense, dust, then we might expect a correlation between the 1.1 mm flux density and the 
total i7ersc/ie/-derived dust column density at that point, or more likely, the product of 
the i/ersc/ie/-derived column density and temperature. This correlation would likely exist 
whether the dust is heated internally by a compact pre-main-sequence object or externally 
by the interstellar radiation field. To test this idea we have plotted in Figure 11 the 1.1 
mm flux in an 80" aperture versus the product of the ifersc/ieZ-derived column density and 
temperature at the position of the 1.1 mm source, averaged over 80". We have not included 
the galaxy 3C111 nor the hot, luminous source associated with LkHalOl whose core dust 
temperature is unlikely to be well-sampled with the Herschel beams and whose flux levels are 



saturated at several Herschel wavelengths. Figure 11 shows a rather good correlation between 
1.1 mm flux and the product of the ifersc/ieZ-derived dust column density and temperature 
for all the other 1.1 mm sources. Thus, the three 1.1 mm sources without associated compact 
70 fim sources probably deserve future investigation as possible pre-stellar cores^j 



4. Individual Sources 

There are several regions within our Herschel maps that are notable for either the 
number of young objects or strong emission at the longer wavelengths. H. Broekhoven-Fiene 
et al. (2013, in preparation) note the large number of Spitzer YSO candidates near LkHalOl 



and in the filament extending north and west of it by roughly 1 degree. Figure [12] shows a 
3-color composite image of this area with the positions of the Herschel sources from Table 
2] marked. Roughly 60% of the likely young far-infrared sources listed in Tableware within 
a 1.2° x 0.5° area (9.4 pc x 3.9 pc) centered on the high-column-density filament shown in 
this area in Figure pi The YSO population in the roughly 4' x 4' core of this region (0.5 



pc square) centered on LkHalOl has been summarized by Herbig et al. (2004) and Andrews 
& Wolk ( |2008 ). More recently Gutermuth et al. (2009) discussed Spitzer observations of 



this cluster in comparison with a number of other young clusters and Wolk et al. (2010) 



have added x-ray data to further define the cluster properties. Our Herschel observations 
are complementary to these studies in that the central few arcminutes of our images are 



I'Bolocam sources 3 and 4 are positionally associated with the bright F star, HD 27214, but the Hipparcos 
parallax for this star of 11.51 milli-arcsec implies that HD 27214 is much closer than the AMC. 
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Fig. 11. — Plot of 1.1 mm flux density in an 80" aperture from Table [3] versus product 
of Herschel dust temperature and column density averaged over the same size aperture. 
LkHalOl and 3C111 are not included. The correlation coefficient for this distribution is 
0.67. The formal uncertainties in the product of dust temperature and column density are 
insignificant in comparison to the uncertainties due to the particular model used. 
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dominated by the diffuse dust emission, presumably heated by the central bright star as 
well as the dense population of lower luminosity stars surrounding it. Beyond a radius of 
~ 3' from LkHalOl, though, we are sensitive to compact thermal emission from individual 
members of the extended YSO population. 

At the southern end of the filament containing LkHalOl is a group of sources that 
includes Source 49, mentioned earlier as one of the objects with a very low T 6o ; in spite of 
being detected easily by Spitzer 's IRAC in all four bands with relatively blue colors. It also 
is extended in the IRAC images on a scale of ~2-3" in the north-south direction. These 
characteristics are consistent with it being a disk viewed edge-on, where the IRAC emission 
arises from scattered light above and below the plane of the disk along low extinction lines 
of sight to the central star; other interpretations are, of course, also possible. An expanded 
view of this source along with six other objects is shown in Figure [i~3|i. The derived column 
density in this area peaks at 4 x 10 22 cm -2 essentially at the position of Source 49, and this 
area is associated with Bolocam Source 18. The SEDs of Source 49 and several other of the 



coldest sources discussed below are shown in Figure [14} As mentioned earlier, for these six 
examples of the coldest source SEDs we have extracted approximate PSF-fit fluxes at 250/im 
which show that the SEDs of all of these objects peak shortward of 250/im. Therefore, we 
have not found any compact objects with extremely cold SED's, although this may be related 
to our requirement for detectable 70/im emission. 

Further north along the filament is a tight collection of young objects that includes one 
of our two 70 (Ltm sources that was undetected by both Spitzer and WISE in the mid-infrared, 
Source 42 in Table [2] with spectral slope a ~ 3 and also associated with Bolocam 1.1 mm 
Source 14. A closeup of this region is shown in Figure [13]d. The image also includes 70 /im 
Sources 35, 38, 39, 40, 41, 43, and 44. 

More than a degree northwest of the filament containing LkHalOl is another region of 
high column density containing four objects from Table [2j Sources 14 - 17. Source 16 is a 
pointlike Class I/O object located ~ 16" west of Source 17 (Class II), the latter of which is 
elongated in the direction of Source 16. Both are also associated with Bolocam Source 7, 
the third brightest 1.1 mm emission region. These two objects are shown in the lower left of 



Figure 13:. In the upper right are Source 15, another very cold object with spectral slope 



a ~ 3, and Source 14, a faint Class I object. 



Finally, in Figure |T3} i we show a tight collection of 7 Class and I sources at the 
northwest end of our mapped area. These sources include an extended object at Herschel 
wavelengths that has been extracted as five separate condensations at 70 /im (Sources 3- 
7), three of which are well isolated Spitzer YSOs from H. Broekhoven-Fiene et al. (2013, 
in preparation). At 160 /mi, however, the five sources are blended into a single elongated 
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Fig. 12. — False color image with 70 /im (blue), 160 /im (green), and 250 fim (red) of the 
LkHctlOl area and star-forming filament to its north. As for Figure [7], the locations of the 
sources in Table [2] are marked with a red 'X', and those that are not in the Spitzer YSO 
list of H. Broekhoven-Fiene et al. (2013, in preparation) are also surrounded with a black 
square. LkHalOl is the brightest object in the image. The second black square down from 
the northern boundary of the image is the location of Source 42, shown in more detail in 



Figure 13 The black square at the bottom of the image marks Source 49, also shown in 



more detail in Figure 13 
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Fig. 13. — False color images of the four areas discussed in ^4] with 24 //m (blue), 70 /zm 
(green), and 160 /im (red). The sources are labeled with their reference numbers from Table 
M (for clarity, sources 4 and 6 are not marked in Panel d) . Column density contours from the 
map shown in Figure [8] are overlaid at values of 0.5, 1.0, 2.0, and 3.0 x 10 22 cm" 2 (N112). 
Panels a-d respectively have angular sizes of 6.4', 5.9', 12.8', and 6.9'. 
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Fig. 14. — Spectral energy distributions of the six coldest sources discussed in £|4} Note, 
sources 4 and 42 have only flux upper limits shortward of 70 fim. 
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structure that peaks on the position of the brightest source at 24 /zm, 70 /im, and 160 /im. 
Sources 2 and 8 from Table [2] are also included in this figure. 



Cloud Mass and Structure 



We have performed a preliminary analysis of the diffuse dust emission as described 
earlier by fitting the 160 /zm-500 /im emission with a simple SED that characterizes the dust 
with a temperature and column density. These results shown in Figure [8] reveal a network 
of narrow filaments characterized by column densities of up to a few x 10 22 cm -2 (N112) 
and temperatures that drop to ~ 10 K from the typical value of order 14-15K in the low- 
density parts of our maps. Many of these filaments are associated with Lynds dark clouds 



as indicated in Figure 1 of Lada, Lombardi & Alves (2009). In an initial effort to quantify 



the differences and similarities between the AMC and other star-forming regions we have 
calculated two quantities discussed by other authors for similar clouds. These quantities 



are the cumulative mass fraction as a function of extinction as already mentioned by Lada, 



Lombardi & Alves (2009) for the AMC and the probability density function for the column 



density as discussed by Schneider et al. (2012) and others 



Figure 15 shows the cumulative mass fraction versus K magnitude extinction A# using 
the same conversion factors as Lada, Lombardi & Alves (2009), 2 x N^ 2 /A^- = 1.67 x 10 22 



cm -2 mag -1 . The figure also shows the Lada et al. distribution as read from their Figure 4 
and a version from our data after smoothing the i/ersc/je/-derived column density with an 80" 
HPW Gaussian to try to duplicate the Lada et al. result. Both our native- resolution function 
as well as the smoothed version are well above the Lada et al. function at all values of Ax- It 
is possible that part of this difference is due to the fact that our observations cover only about 



half of the area mapped by Lada, Lombardi & Alves (2009) that has an extinction above A^- 



= 0.2 mag, though it is difficult to imagine quantitatively how such an areal difference could 
have such a large effect on the derived mass function. Another significant difference, besides 
just the basic technique, is the higher angular resolution of our data compared to the NICER 
optical extinction method, ~ 35" versus ~ 80", but our smoothed mass function is also well 
above that measured with the NICER technique. Yet another possible explanation for the 
difference is that there may exist a population of dust even colder than that sampled by 
the Herschel observations that could be contributing to the extinction-derived dust masses. 
The fact, however, that we have clearly sampled dust down to T~ 10 K and that the diffuse 
dust emission is significantly warmer at temperatures of order 14-15 K argues against this 
hypothesis. We have also investigated whether the difference might be due to differences in 



assumed dust properties. The underlying relationship which would affect Figure 15 is the 
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ratio of far-ir optical depth used to derive our gas column density to the near-ir dust optical 



depth measured by Lada, Lombardi & Alves (2009). If, for example, the K-band extinction 



were smaller at any given value of far-ir optical depth, then the left side (low Ak) part of 
the i/ersc/ie/-derived mass function would shift closer to that found by |Lada, Lombardi"^ 



Alves (2009). On the other hand, however, the right side of our mass function (high Ak) 
would then drop far below the corresponding part of the NICER-derived mass function, 
particularly for the smoothed (80") version of our mass function. Finally, it is also possible 
that the 2MASS-derived NICER extinctions do not sample well the highest column-density 
parts of the cloud. In any case, it will be very interesting to compare these results with 
i/ersc/ie/-derived column density maps for the OMC to see if a similar discrepancy exists in 
that cloud between the NICER and SED-fitting methods. 

The total mass in our mapped area is ~ 4.9 xlO 4 M with 4.89 xlO M above an 
extinction of Ak = 0.1 mag with the various assumptions mentioned earlier. As shown in 
Table |4j these mass values are about a factor of two below those found by Lada, Lombardi "EE 



Alves (2009) who observed a significantly larger area, most of which is occupied by relatively 



low column density material. So within the uncertainties the total masses are in reasonable 
agreement. Table [4] also illustrates numerically the difference in mass distribution in that our 
estimated mass above Ak = 1 mag is three times larger than that found by |Lada, Lombardi 



& Alves (2009) despite the smaller area covered in our study. 



Figure 16 shows the probability density function (PDF) of column density for the AMC. 



This figure also shows three possible fits to portions of the PDF, two log-normal distributions 
for the central part of the PDF, and a power-law falloff for the high-extinction end for 
comparison with other recent studies of Gould Belt clouds. This observed distribution is 
qualitatively similar to other published column density PDFs, but peaks at a moderately 



lower column density, 2x 10 cm , than that for the Aquila region, 4 x 10 cm (Andre 



et al. 2011). The power-law slope of -3 at high extinctions is comparable to that found for 



Aquila, as well as for the Rosette Nebula by Schneider et al. (2012). 



Qualitatively Figure 16 supports the suggestion by Lada, Lombardi & Alves (2009) that 



the AMC has relatively less high-column-density material than more prolific star-forming 
regions like the Orion Molecular Cloud. This conclusion will be able to be further quantified 
as the Herschel data on the OMC become available for comparison, since it is clear already 
from our results that the resolution of the observations and technique used may be important 
to the detailed results. Further analysis will hopefully also provide some insight into the 
underlying physical mechanisms that lead to these differences. 
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Fig. 15. — Normalized cumulative mass fraction for the area observed versus K magnitude 
extinction using the same assumptions as those of Lada, Lombardi & Alves (2009) (solid 
line). The dash-dot line shows the function derived by Lada, Lombardi & Alves (2009) with 
the NICER technique as read from their Figure 4, which is reported as smoothed with an 
80" HPW Gaussian. The dashed line shows the result of smoothing our i/ersc/ie/-derived 
column density map with an 80" HPW Gaussian. 
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Fig. 16. — Normalized probability density function of area versus column density. The 
dashed line and dash-dot line show two different log-normal distributions that can fit portions 
of the PDF. The dotted line shows a power-law slope of -3.1 that roughly fits the extended 
tail of the distribution at high column densities. 
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5.1. Star Formation Versus Column Density 



We have already discussed the fact that virtually all the young stellar objects are found 
along the high column density filamentary structure shown in Figure [8] In fact, if we confine 
our sample to the Class through II SED objects that are likely to be young enough that 
they are still close to their birthplaces in the cloud, there is only one YSO outside the regions 
of the cloud with N#2 < 5 x 10 21 cm" 2 as mentioned earlier. We can investigate whether 
there is a quantitative as well as qualitative correlation by smoothing both distributions and 



comparing the YSO surface density with the gas surface density. Figure [T7| shows the result 
of this comparison, where we plot the surface density of the 68 Class O-II YSOs and the 
i/ersc/ie/-derived gas column density, both smoothed with a 0.2°(1.6 pc) half-power-width 
Gaussian. For the YSOs we used the union of Spitzer (H. Broekhoven-Fiene et al. 2013, 
in preparation) and Herschel (this paper) objects. The two highest concentrations of YSOs 
are found in the LkHctlOl cluster and in a clump in the northern filament about 3/4° north 
of LkHalOl. The derived column density is also highest in these two circcts cts smoothed to 
a 0.2° half-power-width. There is, however, a less perfect correlation between YSO surface 
density and derived gas column density at the intermediate levels, but it is certainly true that 
the greatest number of YSOs are found in the regions with the highest concentration of dust 
and gas. Conversely, in lower column density areas, but still above the general background 
level, essentially no YSOs are found. 



Kennicutt & Evans (2012) have reviewed this subject extensively in the context of both 



Galactic and extragalactic star formation. Likewise Lada, Lombardi & Alves (2010) and 



Heiderman et al. (2010) have attempted to compare star formation rates and gas surface 



densities in multiple local star-forming environments. All of these studies find a roughly 
power-law relation between gas density and star formation rate over some range of gas 
densities. We can quantify our own conclusions by computing the ratio of the two maps 



plotted in Figure [T7| as a function of the derived column density. Figure [P8] shows the average 
of this ratio in ten column density bins expressed as Ax (mag) from the map smoothed with a 
0.2° Gaussian. This plot shows clearly that there is a strong power-law relationship between 
star formation and column density in the AMC. The surface density of young stellar objects 
increases rapidly at the highest column densities. The slope of this relationship is 4.0. This 
conclusion does not depend strongly on the details of our sample since the distribution of 
Spitzer and Herschel YSOs is similar within the small-number statistics, and our derived 



column density map is qualitatively similar to that of Dobashi et al. (2005). This slope 



of 4.0 is comparable to, but slightly less than the slope of 4.6 derived by Heiderman et al. 



(2010) for an ensemble of star-forming regions at comparable levels of extinction (surface 
density). We do not, however, see any obvious break point in the relation between YSO 
surface density and gas column density, though such a break point might be masked by the 
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smoothing process necessary to deal with relatively small number of YSOs. Note also that 
our maximum smoothed extinction is ~ 1.0 mag, equivalent to a visual extinction of 
Ay ~ 8 mag which is just below the level where Heiderman et al. (2010) see a break in their 
power-law slope. 




Fig. 17. — Derived Njj2 column density smoothed with a 0.2° (HPW) Gaussian shown in 
colored contours versus the YSO surface density smoothed in the same way in black dash-dot 
contours. The YSO contours are in steps of 1, 2, 4, and 8 YSOs per smoothed beam for the 
Class II objects. The smoothed column density ranges from 1 to 8.6 xlO 2 
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5.2. Ionizing Environment 

Sharp edges suggesting shaping by photoionization are seen along the southern border 
of the cloud (Fig. [7]), particularly in the west near £ = 160.5°, b = —9.5°, even down to 
low-column densities. The effects of photoionization, however, are modest because there is 
no evidence for a temperature gradient indicative of dust heating. The sharp cloud edges are 
reminiscent of the similar but stronger effects seen in the Oph North region (upper Scorpius) 



which is being photoionized and shaped by the runaway O star ( Ophiuchi (Hatchell et al. 



2012). 



A possible source of photoionization is the 07.5 III star £ Per (£ = 160.4°, b = —13.1°). 
The star £ Per illuminates the California Nebula, a bright infrared nebula located between 



the star and the AMC in projection. Hoogerwerf et al. (2001) model £ Per as a runaway 



O star that was ejected from the Per OB2 association, and that now has a distance of 360 
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Fig. 18. — Plot of the ratio of smoothed YSO surface density to smoothed column density 
(normalized to the peak) relative to the smoothed column density (expressed in K magnitude 
extinction Ak), both from Figure 17 



pc. The proper motion and distance uncertainties, however, are also consistent with the 
interpretation that £ Per is closer to the AMC cloud, and hence able to influence the AMC 
cloud boundary. In any case, the influence of £ Per on the structure of the AMC is minimal, 
and there is no indication of enhanced star formation along the southwest cloud boundary. 



6. Comparison of the AMC and OMC 

As mentioned above, a comparable Herschel study of the OMC is not yet complete so 
it is not possible yet to make a detailed comparison of the star formation and interstellar 
medium between the two giant molecular clouds in the far-infrared. We can, however, discuss 
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briefly the differences in star formation rate on the basis of our deep Herschel observations. 
H. Broekhoven-Fiene et al. (2013, in preparation) have shown that on the basis of Spitzer 
searches for infrared-excess objects, the AMC appears to have about 5% of the YSO popula- 
tion over a similar that of the OMC QMegeath et"aL|2012fl . Our Herschel photometry 
has discovered a few additional young objects on the basis of 70 /im and 160 /xm fluxes, but 
there is clearly not a significant population of deeply buried YSOs. This ratio of 20:1 for 
the star formation rates is, though, not as great as the ratio of the incidence of very massive 
stars. For example, the seminal study of Blaauw (1964) found more than 50 O and early B 
stars in the Orion OB associations, whereas the AMC probably has only one early B star. 
This suggests that whatever the reasons for the lower star formation rate in the AMC despite 
its total mass, the rate for high mass stars is depressed even more in the AMC relative to 
that in the OMC. 



7. Summary 



We have completed the census of star formation in the AMC that began with the Spitzer 
Gould Belt survey of H. Broekhoven-Fiene et al. (2013, in preparation). We have found a 
modest number of additional YSOs, 11, several of which exhibit quite cold SEDs that peak at 
150 /im-200/im. We also mapped a subsection of the Herschel area with Bolocam at 1.1 mm 
and found 18 cold dust sources whose fluxes are well correlated with the dust temperature 
and column density derived at shorter wavelengths with Herschel. We have analyzed the 
distribution of column density and found a strong non-linear relation between column density 
and YSO surface density. We have compared our derived cumulative mass fraction with that 



found by Lada, Lombardi & Alves (2009 ) with the NICER method and noted some differences 



that may be due to a combination of factors including: area covered, angular resolution, and 
details of the methods. The cumulative mass fraction and the probability density function 
for the column density are both qualitatively similar to other clouds for which they have 
been derived but may suggest that the AMC is dominated by lower column density material 



than other clouds with higher rates of star formation as suggested by Lada, Lombardi fc 



Alves (2009). The star formation rate in the AMC appears to be a factor of 20 below that 



in the OMC for typical stars, but an even greater difference exists at the high-mass end of 
the IMF. 
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Table 1: AOR List 



AOR Name 



ObsID 



Field Center 



Comments 



SPParallel-aurwest-orth 1342239276 

SPParallel-aurwest-norm 1342239277 

SPParallel-aurcntr-orth 1342239278 

SPParallel-aurcntr-norm 1342239279 

PPhoto-secluster-orth 1342239441 

PPhoto-secluster-norm 1342239442 

SPParallel-aureast-orth 1342240279 

SPParallel-aureast-norm 1342240314 
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Table 4: Cumulative Mass Versus Extinction 
A K (mag) Mass (this study a ) Mass (Lada09 b ) 



Mag M M 

OTO 4.9 xlO 4 NA 

0.1 4.89 xlO 4 1.12 xlO 5 

0.2 4.28 xlO 4 5.34 xlO 4 

1.0 3.29 xlO 3 1.09 xlO 3 



"Total survey area is 16.5 square degrees. 
b Total survey area is ~80 square degrees. 
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